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We describe the properties of a pair of ultracold bosonic atoms in a one-dimensional harmonic 
trapping potential with a tunable zero-ranged barrier at the trap centre. The full characterisation 
of the ground state is done by calculating the reduced single-particle density, the momentum dis- 
tribution and the two-particle entanglement. We derive several analytical expressions in the limit 
of infinite repulsion (Tonks-Girardeau limit) and extend the treatment to finite interparticle inter- 
actions by numerical solution. As pair interactions in double wells form a fundamental building 
block for many-body systems in periodic potentials, our results have implications for a wide range 
of problems. 
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I. INTRODUCTION 

The last two decades have seen considerable exper- 
imental advancement in the area of cooling and trap- 
ping neutral atoms |1| , with one of the crowning achieve- 
ments being the realization of Bose-Einstein condensa- 
tion (BEC) 2]. There continues, today, intense exper- 
imental investigation into systems of trapped, ultracold 
atoms, with potential deployment of this technology in 
the fields of precision interferometry and quantum in- 
formation processing 3]. In particular, exciting advance- 
ments have been reported on the behaviour of cold atoms 
in periodic potentials, which can be created from the op- 
tical dipole forces arising from several crossed, interfer- 
ing laser beams [J, H[ . Such arrangements have allowed 
experimentalists to trap and control small numbers of 
particles on tightly-confined, individual lattice sites and 
thereby severely restrict their centre-of-mass dynamics. 
Since such potentials can be applied in selective direc- 
tions in space, these techniques allow the creation of ef- 
fectively lower-dimensional systems [(| 0, H, Q • 

A further external handle for control over cold atomic 
many particle systems is the ability to change the inter- 
particle scattering length, allowing access to ideal, as well 
as strongly correlated regimes. This can be accomplished 
by using Feshbach resonances [l(| or by tuning of the ef- 
fective mass of particles moving in a periodic potential 
i- 

Combining these techniques has permitted the experi- 
mental realization of atomic gases in the so-called Tonks- 
Girardeau (TG) regime, wherein a quasi-lD quantum gas 
of strongly-interacting bosons acquires fermionic prop- 
erties [HI, [H, [HI]. Not surprisingly, these experimen- 



tal advancements have motivated much theoretical in- 
vestigation of systems of strongly interacting bosonic 
ID, subject to different confining potentials 



'Electronic address: 
t Electronic address: 



j.f.mccann@qub.ac.uk 



thbusch@phys.ucc.ie 



The theoretical description of a sufficiently dilute sys- 
tem can be achieved by restricting consideration to one- 
and two-particle effects, only. As such the fundamen- 
tal building block for the description of the many-body 
system is the system of two interacting particles, sub- 
ject to some trapping potential. In addition, for the 
low momenta associated with ultracold particles it be- 
comes possible to represent the particle-particle interac- 
tions through a pseudopotential, whereby the descrip- 
tion of the particle-particle interactions depends only on 
the s-wave scattering length [l8[ ■ Previous work has re- 
ported the analytical solution for a pair of particles in 
isotropic (l9l | and anisotropic [2fj| | three-dimensional har- 
monic confining potentials, within the pseudopotential 
approximation. The analytic solution for the ID case is 
also presented in [l9[ . It is straightforward to adapt this 
solution to the problem of a single particle in a (5-split 
harmonic trap t l€3j . The delta-split trap potential may 
be viewed as a generic model for double well situations 
or, alternatively, as a good approximation to the problem 
of a trap with an impurity at the centre. In (l6l | analytic 
ID single-particle eigenstates are used to construct the 
many-body ground state for a system of N particles con- 
fined by a (5-split trap in the TG limit. 

As the single-particle eigenstates are known for arbi- 
trary barrier strength, it is straightforward to obtain an 
analytic expression for the two-particle ground state in 
the TG limit, while for finite interactions a numerical 
scheme is required. In this work we analyse the physics of 
a boson pair including the reduced single-particle density, 
the momentum distributions and the two-particle entan- 
glement, which we quantify by means of the von Neu- 
mann entropy 0, Hi, [IS 0, [S US H3, HI • In particu- 
lar we consider how these properties of the ground state 
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may be altered as both the barrier strength and interac- 
tion strength are varied. Ultra-cold few-boson systems 
in a double-well trap have recently received a thorough 
numerical investigation. In (29l . l30l l3l| the authors em- 
ploy narrow-width Gaussians to model both the contact 
potential and central splitting potential. They proceed 
to use a Hartree-Fock type method to investigate some of 
the many-body properties. The results of our work agree 
well with this numerical approach in limiting cases. 

The remainder of this paper is organized as follows. 
In Sec. |TI] we describe the effective Hamiltonian and the 
assumptions of the model. In Sec. IIIII the single-particle 
eigenstates are reviewed and used in Sec. HVI to construct 
an analytical representation for the two-particle ground 
state in the TG regime (i.e. the Tonks molecule). Us- 
ing this analytical representation we investigate the de- 
pendence on the barrier strength of the reduced single- 
particle density, the momentum distribution and the von 
Neumann entropy for a boson pair. Sec. [V] employs a 
discretization scheme to allow for the variation of the 
interaction strength between the particles. The compu- 
tational method is outlined and a set of results are pre- 
sented. Finally, in Sec. IVII we make some concluding 
remarks and comment on the experimental realization of 
the proposed system. 



II. MODEL HAMILTONIAN 

Consider a system of two identical bosonic atoms which 
are confined in a highly anisotropic harmonic trapping 
potential where the trapping frequency in the perpen- 
dicular directions, is much larger than in the axial 
direction, uj±_ u>. The associated length scales are 
d± = yjh/muj^ and d — ^/i/raw. As a result of the 
large energy level separation, associated with the trans- 
verse eigenstates (Hlu±), at low temperatures the trans- 
verse motion is restricted to the lowest mode. In this 
case the system can be treated as quasi- ID and may be 
described using the effective Hamiltonian, 



H = ^ hi + 9\d5{x 2 - x\) 



(1) 



i=l,2 



where the single-particle Hamiltonian, hi is given by 

h 2 d 2 l 



2m dxj 
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n5(xi) 
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Here, m is the particle mass, and Xi (i = 1,2) is the 
ID position coordinate of particle i. The last term of 
the single-particle Hamiltonian represents a point-like 
barrier located at the origin and the parameter k > 
determines the strength of this barrier. The quan- 
tity gio represents the interaction strength, and is re- 
lated to the ID s-wave scattering length (clid) through 
cjid = — 2h 2 /raaio . In turn, am is related to the actual 
three-dimensional s-wave scattering length, a^o, through 



aio = —d\j2a^u{X ~ Cazp/di ) , where C is a constant 
of value C= 1.4603... $2}. 

In the limit of tight confinement, the free-space pseu- 
dopotential approximation for the par ticle-particle inter- 
actions becomes compromised |33l l34j . In this case, one 
may obtain the eigenenergies for the system by employ- 
ing an energy-dependent scattering length and solving 
for the energy eigenvalues self-consistently [35j, |3g, |37j. 
For current purposes it is supposed that we are in the 
regime for which the pseudopotential approximation is 
still valid and the ID collisional coupling, g\D, acts as 
a parameter for the system. This regime requires that 
the range of the interparticle interaction be much smaller 
than the characteristic length scale of the confining po- 
tential [H, H HI, US HI (i.e. ai D < d). 



III. SINGLE-PARTICLE EIGENSTATES 

We can rewrite the single particle Hamiltonian in 
eq. ([2]) according to the rescaling x = dx, where d is the 
ground state extent in the axial direction, as introduced 
above, 



2dx 2 



— + ~x 2 + kS(x) , 



(3) 



where the scaled barrier strength is now given by k 



(hujdy 



The time-independent Schrodinger equation 



for this system then reads 



h4> n (x) = E n cj> n (x) . 



(4) 



Due to the scaling, the energies E n are given in units of 
fko. At this point, for convenience, we drop the 'bar' on 
all quantities, and acknowledge that we are, henceforth, 
dealing in the scaled quantities just described. The ana- 
lytic solution to eq. ([4]), for those eigenfunctions of even 
symmetry, can be found as, (l9| . 



-u 



~2~ 



0,2,4. 



(5) 

Here Af n is the normalization constant and U(a, b, z) are 
the Kummer functions [38j |. The corresponding eigenen- 
ergies, En, are determined by the roots of the implicit 
relation, fl9( ] 



- k= 2 
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(6) 



By contrast, the antisymmetric eigenfunctions vanish at 
the origin and are unaffected by the barrier. They are 
therefore given by the odd eigenstates of the unperturbed 
harmonic potential (k = 0) 



4> n (x) = N n H n (x)e x 2 n= 1,3,5. 



(7) 



where H n (x) is the n th order Hermite polynomial. The 
corresponding energies are given by the eigenvalues of 



3 



(a) k = 




(b)K=l 




(d) K — 00 




the odd parity states of the harmonic oscillator, E n 
2) 



n - 



Considering eq. (J6J) in the limit k — > 0, we find 
E n = i, §,§,..., and the even eigenstates are simply 
given by the even harmonic oscillator solutions. On the 
other hand, for k — > 00, these energies converge towards 
E n = |, |, tt, •• ., and each even eigenstate becomes de- 
generate with the next highest-lying odd parity state. 



IV. TONKS MOLECULE 

In the limit gm — > 00 the point-like, impenetrable, 
interaction between the two atoms can be represented as 
a constraint on the allowed bosonic wavefunction, \ff 

El BE [3, 



k I 



tyf(xi,x 2 ) = if xi = £2 for all fe 



(8) 



where A; is an index labelling the eigenstates. One can see 
immediately that this constraint is equivalent to the ex- 
clusion principle for a corresponding system of two spin- 
aligned fermions, which is a symmetry that gives rise to 
the Bose-Fermi mapping theorem [12, HH, G3 ■ It allows 
one to solve the strongly interacting system of two bosons 
by solving the, often more accessible, system of two non- 
interacting fermions, then properly symmetrising the fi- 
nal wavefunction. In particular, the ground state of the 
two-boson system (the Tonks molecule), , is related 
to the non-interacting fermionic ground state, , by 



^(x u x 2 ) = |*o 0*1 > x *)\ 



(9) 



The fermionic ground state, (xi,x 2 ), is given by the 
Slater determinant of the two lowest single-particle or- 
bitals, so that, 
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\<t>o(x\)(t>i(x 2 ) - 0o(x2)0i(ari)| 



= N_ e -(x?+x?) 
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where Af is the normalization factor. 

Fig. Q] shows the two-particle wavefunction for the 
Tonks molecule in (£1,2:2) space, given by eq. p0[) . for 
different values of k. For k = (Fig. Ufa)) the nodal line 
along xi = x 2 reflects the infinite repulsion of the TG 
limit, or equivalently the exclusion principle of eq. (|8j). 
The distribution of the two-particle wavefunction shows 
a strong correlation between an x\ > and an x 2 < 
coordinate, and vice versa. Increasing ft to 1, 2 and 
00 (Figs. QJb),(c) and (d)), the wavefunction is reduced 
along the lines xi — and x 2 — due to the strength- 
ening potential barrier at the origin. In this process the 
wavefunction also becomes increasingly squeezed along 
the line xi — —x 2 , indicating the localization of one par- 
ticle on each side of the barrier. We note that for values 



FIG. 1: (Color online) Ground state wavefunction for a boson 
pair in a harmonic trap with a 5-barrier along xi = X2 = 0, 
of strength (a) k = 0, (b) ft = 1, (c) k = 2 and (d) n — 00. 
The corresponding, scaled ground state energies are Eo = 
2.0, 2.4, 2.6 and 3.0, respectively. In each plot the horizontal 
and vertical axes run from -6 to +6, in scaled units. 



of k > 2 there is no appreciable change in the two-particle 
density with barrier strength. 



A. Reduced single-particle density 

A quantity of fundamental importance in many-body 
physics is the reduced single-particle density (RSPD), 
given by, 



p{x,x') 



^{x,x 2 )^{x' ,x 2 )dx 2 



(11) 



The RSPD is illustrated in Fig. [U for four different 
barrier strengths, corresponding to the same values ex- 
amined in Fig. [U i.e. k — (a), n — 1 (b), K = 2 (c) and 
k = 00 (d) . The RSPD expresses the self correlation and 
one can view p(x,x') as the probability that, having de- 
tected the particle at position x, a second measurement, 
immediately following the first, will find the particle at 
the point x' . Classically, p(x,x') = 5(x — a;'), and one 
can see from Fig. [5]that a strong enhancement of p(x, x') 
exists along the line x — x' . In the absence of any bar- 
rier, Fig. Ella), the significant off-diagonal contributions 
reflect the delocalization of an individual particle, since 
there is a non- vanishing probability that the second mea- 
surement, x' , may find the particle anywhere in the trap. 
Increasing the barrier strength, as seen in Figs.[T]Jb), (c) 
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(b)K=l 




FIG. 2: (Color online) Reduced single-particle density matrix 
for the Tonks molecule, for barrier strength (a) k — 0, (b) 
k = 1 , (c) k = 2 and (d) k = oo. Each plot spans the range 
—6 < x, x' < 6. 



and (d), leads to the emergence of a quadrant separa- 
tion. For a stronger barrier the contributions in the off- 
diagonal quadrants diminish. In particular, in Fig. [TJd) 
these off-diagonal contributions to p(x, x 1 ) vanish alto- 
gether. The strong barrier restricts tunnelling from the 
left side of the well to the right, and vice versa. In this 
scenario, the ground state of the system is comprised of 
each member of the boson pair in a separate half-well. 



B. 



Momentum distribution 



While, due to the Bose-Fermi mapping theorem, the 
density distributions of a sample of bosons and fermions 
becomes identical in the TG limit, the momentum distri- 
bution can still be used for distinction [lj| . The recipro- 
cal momentum distribution, n(k), is calculated from the 
reduced single-particle density 



p{x, x')e- %k ^- x 'Uxdx' , 

OO 

(12) 

Equivalently, one may obtain 
the momentum distribution for this system by consider- 
ing the diagonalization of p(x, x'). The eigenvalue equa- 
tion to be solved is 



i(k) = (2tt)" 



where f_°° n(k)dk = 1 . 



p{x,x')ipi{x')dx' = Xitpi(x) , 



(13) 



where the eigenvalue, Ai, represents the fractional popu- 
lation of the 'natural orbital' ipi(x) such that J^i At = 1. 
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FIG. 3: (Color online) Momentum distribution, n(k), for 
k = 0, 1, 5 and 10, with the corresponding normalised ground 
state energies Eo = 2, 2.4, 2.8 and 2.9. The momentum dis- 
tributions broaden for increased barrier strength due to the 
localization of the particles in separate halves of the trap. 
Also shown is the momentum distribution for the k = oo 
case, given by eq. (|16|l . for which Eo = 3.0. 



Using a discretized form for the quadrature allows one 
to rewrite the integral equation Q13[) as a linear algebraic 
equation. The momentum distribution, n{k), may then 
be obtained from the relation 



i{k) = J £\ i \ fM (k)f 



(14) 



where pi(k) denotes the Fourier transform of the natural 
orbital ipi{x), 



p%{k) 



tpi(x)e lkx dx 



(15) 



Fig. [3] shows the momentum distribution in the TG limit 
for four different values of the barrier strength, k = 0, 1, 5 
and 10. As the barrier strength is increased the momen- 
tum distribution becomes broader. This observation is 
consistent with the earlier observation that the two sep- 
arate particles become individually localized in the two 
separate half-wells. 

In the limit of infinite barrier strength (k = oo) the 
system becomes doubly degenerate which allows us to 
calculate an analytical expression for the momentum dis- 
tribution. By changing the computational basis and 
defining rj(x) = [<t)o{ x ) + ^lO^L where (f> and 4>i are 
the ground and first-excited cigenfunctions of the single- 
particle Hamiltonian ([3]), with 4>q(x) = \(f>i(x)\, we find 
that the wavefunction is only finite in the region x > 0. 
The momentum distribution is the given by the direct 
Fourier transform of rj(x) 



n (k) = — 

7T2 



1 



k 2 e- 



1 3 V 

2'2'2 



1 2 



+ 2 k 6 



(16) 
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This analytic momentum distribution, for the case n = 
oo, is also plotted in Fig. [3] It can be seen that, in the 
limit of large k, the momentum distribution calculated 
from the diagonalisation of the RSPD matrix tends to- 
wards the profile given by eq. (fT6|) . 

C. Ground state entropy 

Entanglement is not only a fundamental quantity in 
quantum mechanics, it is also one of the most impor- 
tant resources in quantum information theory, where it 
is often responsible for the increased efficiency of quan- 
tum algorithms over their classical counterparts. Pre- 
vious authors have shown that the von Neumann en- 
tropy is a good measure of entanglement for a system 
of two bosons [H, [24], [25[ . In the case of indistinguish- 
able particles , however, differentiating between entan- 
gled and non-entangled states requires that one considers, 
simultaneously, both the von Neumann entropy of the 
reduced single-particle density and the Schmidt number 
dHHHI. The Schmidt number is given by the number 
of non-zero eigenvalues, Aj, of the reduced single-particle 
density, p (see cq. (O). In this work we shall use the 
von Neumann entropy to quantify the entanglement in 
the position coordinates, x\ and x 2 , of the boson pair 
and the Schmidt number shall only be discussed when it 
affects the interpretation of the results presented. 

The von Neumann entropy, S, is defined by 

S=-J2^og 2 K- (17) 

i 

and we calculate the values for A^ by numerically diag- 
onalising the RSPD matrix as a function of k. The re- 
sults are shown in Fig. 2J Interestingly, one sees that 
the entropy begins at a value of about 0.985 for k = 0, 
which agrees well with the limiting value suggested in 
(2^ . As k increases, S is seen to increase through a 
value of unity. It peaks for n w 3.4 (corresponding to 
a ground state energy of Eq ss 2.85) before dropping off 
and tending towards unity in the limit k — > 00, corre- 
sponding to a non-entangled state. Identification of this 
state as non-entangled follows from the fact that the von 
Neumann entropy for this state (with an infinite barrier) 
equals unity and the Schmidt number is found to equal 
2, [2^, [27], [28[ . In this situation, the ground state of the 
boson pair is comprised of one particle residing in the 
left half-well and one in the right. However, owing to 
the indistinguishability of the particles, one cannot say 
which particle resides to the left and which to the right. 
This lack of information, arising solely from the indistin- 
guishability of the particles, leads to the value of 1 for 
the von Neumann entropy. Pure statistical correlations 
are of little intrinsic value to any quantum information 
protocol, and the state is regarded as non-entangled. 

By contrast, the point at which S — 1 for the finite 
value of k ~ 1.33, represents an entangled state. This is 
due to the fact that the Schmidt number at this point is 
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FIG. 4: (Color online) Von Neumann entropy, S, for the 
Tonks molecule, as a function of the barrier strength, ft. When 
k — then S w 0.985. As the barrier is strengthened the 
entropy increases to a maximum at k w 3.4. In the limit 
of k — ► 00 then S 1, corresponding to a non-entangled 
state. The bar charts show the values of the Schmidt num- 
bers at the point where (« = 1.33, 5 = 1) and for the limit 
(ft = 00, S = 1). 



> 2, allowing one to classify the state as truly entangled, 
beyond purely statistical correlations, [2ol I27I |28| . 

V. VARIABLE INTERACTION STRENGTH 

For finite particle interactions no analytical solution to 
the inhomogeneous two-particle problem is known (ex- 
cept in the case of k = 0, [Hj]). In this section we, there- 
fore, use a numerical discretization scheme to study the 
ground state properties of the boson dimer as a func- 
tion of varying interaction strength, as well as barrier 
strength. 

Discretization of the spatial coordinates X\ and x 2 is 
achieved by means of a discrete variable representation 
(DVR), Ho, El|. The two-particle wavefunction is repre- 
sented by the direct product 

N 

®(x u x 2 ) = *ij7i(*i)/j(a'a) ■ (18) 

Here is the value of the two-particle wavefunction at 
the mesh point (x% = qi, x 2 — qj), with i,j = 1, 2, . . . , N. 
Clearly, these mesh points are finite in number and will 
be restricted to some region in (xi, x 2 ) space, defined by 
the boundaries a and 6, such that 

a<q t <b i = l,2,...,N. (19) 

The values, ^ij, play the role of variational parameters to 
be found and the /, (q) are a set of N Lagrange functions 
which have the property that they are localized about 
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the mesh points q±, q%, . . . , q^. In addition to satisfying 
the usual interpolation conditions, 



g =500 

° D 



V i,j 



(20) 



one also requires that these Lagrange functions satisfy 
the orthogonality condition 



f*(q)fj(q)dq = \h 



(21) 



Here \ are the generalised Christoffel numbers associ- 
ated with the mesh, Eel, and 



Ai = 1 V i , 



(22) 



for the Cartesian mesh considered in this work. For this 
Cartesian mesh the Lagrange functions are given by 



Mo) 



1 sin[7r(<7 — i)] 
N sm[ir(q - i)/N] 



(23) 



Using the basis expansion of eq. Ijl8|> in the Schrodinger 
equation Q results in a discrete eigenvalue problem that 
can be solved using standard linear algebra techniques. 



A. Reduced single- particle density 

We have calculated p(x,x') using N = 81 mesh points 
in each dimension (i.e. x\ and X2) and a mesh spacing 
of Ax = 0.16. Color density plots of the reduced single- 
particle density are presented in Fig. [5] for four different 
values of interaction strength and four different values of 
barrier strength. Each row illustrates the transition from 
a non-interacting pair (giD = 0) to a strongly interact- 
ing dimer (gm = 500), and each column illustrates the 
transition from a single well (k — 0) to an, essentially, 
split trap (k = 10). 

In the first column of Fig. [5] the non-interacting limit 
(91D — 0) is considered. The increased barrier strength 
at the origin manifests itself by diminishing p(x, x') along 
the lines x = and x' = 0, thus partitioning the struc- 
ture into four quadrants. The even division of p(x, x') 
over all four quadrants reflects the derealization of each 
individual particle over the two half-wells. 

The second column of Fig. [5] shows the same color 
density plots for k — 0,1,2 and 10 for a finite interac- 
tions strength of g\o = 1. In the absence of a barrier 
(k = 0) the RSPD exhibits similar features to the non- 
interacting case, although it expands slightly in both x 
and x' . Strengthening the barrier again gives rise to a 
quadrant structure. However, the presence of repulsion 
reduces p(x, x') in the off-diagonal quadrants, meaning 
that the initial detection of a particle in the left half-well 
precludes its subsequent detection in the right half-well 
and vice versa. Analogously to the Bose-Hubbard model, 
the system will be governed by the interplay between the 
tunnelling (determined by the strength of the barrier, k) 



K=2 



K= 10 




FIG. 5: (Color online) RSPD, p(x,x') as a function of in- 
teraction strength gm =0,1,5 and 500 and barrier strength 
k = 0, 1, 2 and 10. Each individual plot spans the range 
-6.4 < x, x' < 6.4. 



and the on-site interaction (determined by the interac- 
tion parameter, gm)- For a strong barrier (k = 10) and 
finite interaction, there is a blockade and the insulator 
state dominates, with one boson in each half-well. In 
terms of the reduced single-particle density, p(x,x'), this 
leads to the vanishing of the off-diagonal contributions 
as tunnelling of a given particle between the two half- 
wells becomes increasingly unlikely. This behaviour is 
increasingly visible in the third and fourth column when 
the interaction strength is increased to gm = 5 and 500, 
respectively. 

As the interaction strength increases, the first plot in 
the third column shows a clear deviation from the circular 
structure observed in the k = case for lower interaction 
strength. The distribution is now clearly enhanced along 
the line x = x' and reduced in the direction orthogonal 
to this. The stronger repulsive interaction has the effect 
of reducing the 'derealization' of the particles. As the 
barrier strength is increased the off-diagonal contribu- 
tions die-off faster than in the case g±o = L This is due 
to the stronger interactions encouraging the localization 
of the particles at even smaller barrier strengths. The 
superfluid character, that is indicated by the quadrant 
structure, decays already for smaller values of,K. As in 
the case of gio — 1, as the barrier strength is, increased 
one observes the reduction in the off-diagonal contribu- 
tions and in the limit of large k one observes the, almost 
perfect, localization of the two particles in the two sepa- 
rate half-wells. 

Finally, the last column illustrates p(x,x') for very 
strong repulsion, gi£> = 500. As one expects, the re- 
duced single-particle densities closely resemble the plots 
displayed in Fig. [2] for the Tonks molecule. 
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FIG. 6: (Color online) Momentum distributions for varying 
interparticle interaction strength gm = (a), 1 (b), 5 (c) 
and 500 (d). Within each plot the distribution is considered 
for different values of the barrier strength: k = (solid line, 
black), 1 (dashed line, red), 5 (dash-dot line, green), 10 (dot- 
ted line, blue). These calculations have been carried out by 
means of the DVR discretization of the spatial coordinates x\ 
and X2, with N = 61 DVR mesh points in each dimension 
and a scale factor of Ax = 0.16. 



B. 



Momentum distribution 



The momentum distributions, n(k), can be obtained 
from the reduced single-particle density, p(x,x'), using 
the same methods outlined in Sec. IIVBI 

Fig. [Ha) shows the distribution obtained for two 
non-interacting particles (giD — 0) for varying k. In 
this case, the momentum distribution is given by the 
square of the single-particle wavefunction in momentum 
space, |x(fc)| 2 . In the limit of an infinitely strong bar- 
rier the single-particle wavefunction becomes (f>o(x) = 
{2/^) 1 / 2 \x\e- x I 2 and we can calculate the momentum 
distribution analytically, 



n(k) 



r 3/2 



1 - k 2 e~-M 



2 2 



(24) 



The interplay between the first and second terms give rise 
to the secondary peaks seen in Fig. [6]^a) for k — 2,5 and 
10. Physically, these peaks arise due to the interference of 
the particle, split between the two separate half-wells, in 
analogy with a double-slit arrangement. When increasing 
the interaction strength (Figs. Ob)-(d)) these secondary 
peaks disappear, which can be attributed to the increased 
localization of the particles. 

At the same time as observing the emergence of 
these secondary peaks with increased barrier strength for 
gio = 0, one also observes a narrowing of the central 
peak. Strengthening of the barrier causes the ground 
state to shift upwards in energy, this shift will be accom- 
panied by a spreading of the single-particle wavefunction 
in position space, which in turn gives rise to a reciprocal 



narrowing in momentum space. 

Figs. [He) and (d) display the momentum distribution 
in the limit of strong repulsive interactions with the same 
basic trends being observed in both plots. Once again, as 
with the reduced single-particle density, this fact suggests 
that the behaviour of the ground state remains fairly con- 
stant for interaction coupling g\r> > 5, such that these 
finite values of interaction coupling will lead to behaviour 
which is qualitatively similar to the regime of infinite 
repulsive interaction. The results presented for strong 
repulsive interaction (gm = 500) are expected to corre- 
late closely with the momentum distribution obtained in 
Sec. lIVBl for the Tonks molecule, and a detailed compar- 
ison of Figs. [3] and [S^d) verify that this is the case. 

In the large interactions limit the momentum distri- 
bution for k = is observably different from the non- 
interacting case (solid, black lines in Figs.[HIa) and (d)). 
In particular non-Gaussian wings extending to higher k- 
value are observed in the TG regime. One may consider 
the trapping potential to be switched off suddenly, and 
the two-particle wavefunction allowed to expand freely. 
In this case the wavefunction in coordinate space will 
map on to that in momentum space, in the far field limit. 
Clearly, the strong repulsion between the particles in the 
TG limit will lead to a proportion of the ensemble mutu- 
ally recoiling at high speeds and in this way accounting 
for these high-fc wings in the momentum distribution. In- 
creasing the strength of the barrier then has the effect of 
broadening the momentum distribution as the individual 
particles become localized to individual sides of the trap. 
The wings in the momentum distribution are, therefore, 
a signal of a transition into a Mott-insulator type state. 
The spatial localization of the particles is accompanied 
by a broadening in the momentum distribution and this 
is the broadening observed in Figs. [He) and (d). 



C. Ground state entropy 

Variation of entropy with interaction strength 

Let us first examine how the entropy of the two-particle 
system varies as one changes the interaction strength be- 
tween the particles. Similar calculations have been car- 
ried out by other authors, 23], though restricted to a 
harmonic trap without a barrier. This case is represented 
as the lowest (black) line in Fig. [7] For gm = entangle- 
ment is absent and S = 0. As the interaction strength is 
increased the entanglement increases. For gio — > oo the 
entropy saturates at a value of S ~ 0.985 in the absence 
of any barrier, [231 ]. Also shown in Fig. [7] is the variation 
in the entropy with interaction strength when one intro- 
duces a 5-barrier at the well centre. Four different barrier 
strengths are plotted: k = 1 (dashed line), 2 (dash-dot 
line), 5 (dash-dot-dot line) and 10 (dotted line). One 
striking behaviour is noted, as one increases the barrier 
strength, the sensitivity of S to small changes in <?id, 
about giu — 0, is dramatically increased. As the barrier 
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FIG. 7: (Color online) Effect of varying interaction strength 
(gm) on the von Neumann entropy (S) of the ground state. 
The strength of the <5-barrier is taken to be k = (thick solid 
line), 1 (dashed line), 2 (dash-dot line), 5 (dash-dot-dot line) 
and 10 (dotted line). It is seen that for increased strength of 
the central barrier, the entropy shows an increased sensitivity 
to the interaction parameter about the value gm = 0. 



K 

FIG. 8: (Color online) Effect of varying barrier strength on 
the von Neumann entropy, S, of the ground state. The 
strength of the inter-particle interaction is set to gm = 1 
(solid line), 2 (dashed line), 5 (dash-dot line) and 500 (dotted 
line). The initial value of the entropy (i.e. in the absence of 
any barrier) is dictated by the strength of the interparticle 
interaction, with larger interaction leading to increased en- 
tropy. One observes that in all cases, in the limit of a strong 
barrier, the von Neumann entropy saturates at a value S = 1. 



strength is increased the harmonic trap is split into two 
half- wells and the tunnelling between these two half- wells 
is made increasingly unfeasible. As a consequence, the 
ground state of the two-particle system is less capable of 
adapting to changes in the interaction strength between 
the particles, leading to an increased sensitivity of the 
entropy in this respect. For all values of k, as gm — > oo 
the entropy tends to a value close to unity (see Fig. 0|. 



Variation of entropy with barrier strength 

Finally, we present the results of how the entropy of 
the two-particle system changes with the strength of the 
central barrier in Fig. [5] For the case of zero interactions 
the entropy remains zero for all barrier strengths. In the 
presence of a finite interaction the entropy begins with a 
non-zero value, representing the entanglement in the har- 
monic trap with no barrier. As the barrier is strength- 
ened the entropy increases gradually towards unity, and 
saturates at this value. As was discussed in Sec. IIVC1 in 
the limit of infinite barrier strength, the two-particle sys- 
tem will become non-entangled. This is due to the fact 
that, in order to minimize the energy of the system, the 
repulsively interacting particles will localize on opposite 
sides of the well. The only correlations that then exist 
between the particles can be attributed to their indistin- 
guishable nature. 

In the limit of vanishing barrier strength, the larger 
the interaction strength, the larger is the initial value 
of the entropy. As a consequence, for larger values of 
g±D the entropy changes less dramatically as the bar- 



rier strength is increased. It is noted that for the case 
of gm — 500 (dotted line) we are effectively consider- 
ing the TG regime. From Fig. [5] it appears that in this 
regime the entropy remains close to unity for all values of 
k. However, closer inspection of these numerical results 
reveals that this curve actually follows the same trend as 
illustrated in Fig. 2J obtained from the analytical treat- 
ment of the Tonks molecule. This further illustrates the 
correspondence of these DVR mesh calculations to the 
analytical TG results in the limit of infinite repulsive in- 
teractions. 



VI. CONCLUSIONS 

In the present work we have carried out a detailed 
examination of the ground state for two particles in a 
(5-split harmonic trap. We have found that in the pres- 
ence of interactions the reduced single-particle density 
exhibits vanishing contributions in the off-diagonal quad- 
rants in the limit of increasing barrier strength. This 
feature is attributed to the localization of individual par- 
ticles on either side of the split trap, a situation analo- 
gous to the Mott insulator regime in lattice studies and 
also reflected in the corresponding momentum distribu- 
tions. More specifically, in the non-interacting case with 
a strong barrier one observes secondary peaks in the mo- 
mentum density, attributed to interference. These sec- 
ondary peaks vanish in the presence of interaction owing 
to the localization of individual particles. In the Tonks- 
Girardeau limit, increasing the barrier strength has the 
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effect of broadening the momentum distribution, a fea- 
ture that may be explained in terms of the squeezing of 
the wavefunction for the system in position space. Fi- 
nally, we have shown that the von Neumann entropy for 
this system is sensitive to the two parameters of interac- 
tion strength and barrier strength. For a given barrier 
strength, an increasingly repulsive interaction strength 
will cause the von Neumann entropy to saturate at a 
value close to unity. It is found that increasing the 
strength of the barrier has the effect of making the von 
Neumann entropy increasingly sensitive to small changes 
in the interaction coupling about the value of zero cou- 
pling. At the same time, for a fixed value of interaction 
strength, increasing the barrier strength has the effect of 
increasing the entropy of the system. In the limit n — > oo 
the entropy saturates at a value of unity. 

We would like to remark that even though our anal- 
ysis makes use of an idealised 5-function potential, such 
an approximation is known to, not only, encapsulate the 



basic physics, but can also be a very good approxima- 
tion to experimental setups. These include wide traps 
that are pierced by a highly-focused laser beam as well 
as, for example, the situation where a single particle of 
a different species is confined in the centre of the trap. 
Due to the low temperatures, the interaction with such 
a 'quantum dot' would be well described by a point-like 
potential. 
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